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We construct the ID Z 3 parafermionic model based on the solution of Yang-Baxter equation and ex¬ 
press the model by three types of fermions. It is shown that the Z 3 parafermionic chain possesses both 
triple degenerate ground states and non-trivlal topological winding number. Hence, the Z 3 parafermionic 
model is a direct generalization of ID Z 2 Kitaev model. Both the Z 2 and Z 3 model can be obtained from 
Yang-Baxter equation. On the other hand, to show the algebra of parafermionic tripling intuitively, we 
define a new 3-body Hamiltonian H 122 , based on Yang-Baxter equation. Different from the Majorana 
doubling, the H 122 , holds triple degeneracy at each of energy levels. The triple degeneracy is protected 
by two symmetry operators of the system, w-parity P(ijj = e'~) and emergent parafermionic operator 
r, which are the generalizations of parity Pm and emergent Majorana operator in Lee-Wliczek model, 
respectively. Both the Z 3 parafermionic model and jfli 23 can be viewed as SU(3) models in color space. 
In comparison with the Majorana models for SU(2), it turns out that the SU(3) models are truly the 
generalization of Majorana models resultant from Yang-Baxter equation. 

PACS numbers: 02.20.Sv, 03.67.Lx, 03.65.Fd, 75.10.Pq 


The double degeneracy of a pair of Majorana zero modes in condensed matter system has attracted much attentions due to its 
potential applications in quantum computation and quantum information [ 1-6]. It is well known that this topologically protected 
doubling is immune to local perturbations. Taking ID p-wave Kitaev model [1] as an example, the Majorana mode appears in 
the topological phase where the two free Majorana fermions 71 and 72Ar can be excited without cost of energy at the two ends 
of the chain model and compose a non-local complex fermion, hence the ground state possesses double degeneracy. The two 
degenerate states can be differentiated by the electron number parity operator Pm = i-e-, one state possesses parity 

— 1 with odd electron occupation number, while the other possesses parity +1 with even electron occupation number. On the 
other hand, to give an intuitive analysis about the Majorana doubling, Lee and Wilczek [7] proposed a 3-body Hamiltonian 
Hm = 1(0:7172+/37273+K7173), where the symmetry operators Pm and emergent Majorana operator Tm lead to the doubling 
at any energy level. 

In our previous paper, we have shown that both the Kitaev model and the Lee-Wilczek model can be derived from the 
4x4 matrix representation of Yang-Baxter equation(YBE) [8]. The applications of YBE [9-14] in constructing many body 
Hamiltonian had been discussed in various papers [15-18]. Specifically, based on the Majorana representation of Yang-Baxter 
equation Ri{9) = e®2'i7i+i ^ ^ve take the time derivative of 9 in Ri{9) to obtain the ID Kitaev model [8]. To self-contain, we hrst 
recall some results related to R\{9) = which emerges from the 4x4 matrix representation of YBE, which reads 


P'(0) = 


cos 6* 0 0 sin 9 ' 

0 cos 9 sin 9 0 

0 — sin0 cos 9 0 

. — sin 9 0 0 cos 9 _ 




When 9 = f, the P'-matrix turns into the braid operator P'[19] 
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( 2 ) 


In our previous paper [8], we have shown that the Majorana representation relates to the 4x4 matrix representation of YBE 
R{{9) = in tensor product space through Jordan-Wigner(J-W) transformation, here ax and ay are Pauli matrices, i 

and i-bl signify lattice sites. The J-W transformation transforms spin-2 operators at lattice sites into spinless fermions through 
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where tr^ are spin ladder operators, and a„ are spinless fermions. Define the Majorana fermion[l] 

72i-i = al+Qj, 72j = i(a]-aj), {7*,7j} = (4) 

Substituting equation (3) and (4) into equation ( 1), we can express equation ( 1) as next nearest neighbor interaction of Majorana 
operators, 

(5) 

where 721-1 satisfy Clifford algebra { 721 - 1 , 72 j-i} = 2i5ij. Based on the Clifford algebra, the Rl{0) can be rewritten as 

( 6 ) 

It is easy to check that the Ri (0) satisfies YBE. Hence the matrix representation of solution of YBE and Majorana representation 
of braid operators are well related. 

A question is raised naturally. One spin site corresponds to 2 subcells of Majorana fermions that are related to 4x4 YBE in 
the tensor space due to the 4-d representation of Temperley-Lieb algebra. On the other hand, the 9x9 form of solution of YBE 
has been known, then could we extend the above discussion to the new type of “Majorana fermions” with 3 subcells on one spin 
site? The answer is yes. Instead of SU(2), the SU(3) operators should naturally be introduced. Eor the convenience, we call the 
space color space. 

Since the Majorana models hold 2-fold degeneracy, for SU(3), how can we extend the Majorana double degeneracy to triple 
degeneracy? In other words, can we construct the extended ID Kitaev model holding triple degenerate ground state? Indeed, 
similar to those in constructing Majorana models via 4x4 matrix solution of YBE, we can find the triple degenerate models 
based on the 9x9 matrix representation of YBE [20]. 

In this paper, we make the following progress: 1) Based on the 9x9 matrix representation of YBE and the 3x3 3-cyclic 
representation of SU(3) generators (see Supplementary), we make the decomposition of the 9x9 matrix by tensor products of 3- 
dimensional matrices. By defining generalized SU(3) J-W transformation, we transform the SU(3) sites into non-local operators 
and obtain the new representation of YBE. 2) We obtain the Z 3 parafermionic chain with triple degenerate ground states in color 
space and express the chain with three types of fermions, besides that, the Zn case is discussed; 3) In Z 3 parafermionic model, 
the topological phase transition is signified by the triple degeneracy of ground states and the topological winding number; 4) To 
give an intuitive explanation of the triple degeneracy and analyse the algebraic structure in it, we construct a 3-body Hamiltonian 
and find its symmetry operators that lead to the tripling. 


IfESULTS 


Review of two Majorana models. To preserve the self-consistency of this paper, firstly, let us give a brief introduction to the 
construction of Majorana models based on YBE. The intrinsic connection between the solution Ri{9) = of YBE and 

Kitaev model is that both of them possess Z 2 symmetry. Next we review the Kitaev model derived from YBE. 

We imagine that a unitary evolution is governed by Ri{0). If only 9 (tan0 is the velocity m of a particle) in unitary op¬ 
erator Ri{9) is time-dependent, we can express a state \il>(t)) as = Ri{9{t))\ij}(0)). Taking the Schrodinger equation 

= H{t)\ip(t)) into account, one obtains: 

m^[i?ilV'(o))] = iT(f)4|^(o)). (7) 

Then the Hamiltonian Hi{t) related to the unitary operator Ri{9) is given by: 

H,{t)=ih^Rr\ ( 8 ) 

Substituting Ri{9) = into equation (8), we have 

Hi{t) = (9) 

If we only consider the nearest-neighbour interactions between two Majorana fermions(MEs) and extend equation (9) to an 
inhomogeneous chain with 2N sites, the derived chain model is expressed as [8]: 

N 

Hy^ — ih^^{9i'y2k-ll2k+92l2kl2k+l), ( 10 ) 

fe=l 

with 9i and 02 describing odd-even and even-odd pairs, respectively. This is exactly the Kitaev model derived from YBE. 

The properties of ID Kitaev model are well known: 
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1. In the case 0i > 0, 02 = 0, the Hamiltonian reads: 

N 

Hiriv = i^^^l72fc-l72fc- (11) 

k 

As defined in equation (4), the Majorana operators 72 /c-i and 72 ^ come from the same ordinary fermion site k, 
i 72 /c-i 72 fc = 2 a 7 fe~l («! and au are spinless ordinary fermion operators). Hi simply means the total occupancy of 
ordinary fermions in the chain and has U(l) symmetry, aj e'^aj. The ground state represents the ordinary fermion 
occupation number 0. This Hamiltonian corresponds to the trivial case of Kitaev’s. 

2. In the case 9i = 0, 02 > 0, the Hamiltonian reads: 

N 

Htopo = i/iy^ 0272fc72/c+l- (12) 

k 

This Hamiltonian corresponds to the topological phase of ID Kitaev model and has Z 2 symmetry, aj —o-j- Here 
the operators 71 and 72 Ar are absent in H 2 . The Hamiltonian has two degenerate ground state, |0) and |1) = dl|0), 
= ( 71 —i 72 tv)/ 2 . This mode is the so-called Majorana mode in ID Kitaev model. 

On the other hand, as pointed out by Lee and Wilczek in Ref. [7], the double degeneracy of Majorana models Hm is due to 
two symmetry operators, the parity operator Pm and emergent Majorana operator Tm. For instance, in 3-body Majorana model 
with the Hamiltonian 


Pm = i (q:7i72-I-/37273+k7i73) , (13) 

the symmetry operators and commutation relations are 

Pm = (-l7^ Tm =-i7i7273; (14) 

[Pm, Pm] = 0, [PmjTm] = 0, {TmjPm} = 0. (15) 

Clearly, in the basis of both Pm and Pm are diagonal, Tm transforms the states with Pm = ±1 into the states with Pm = t 1. 
Therefore the Hamiltonian possesses Majorana doubling. 


Yang-Baxter equation and 3-cyclic SU(3) generators. Since the YBE is properly applied in constructing Z 2 Kitaev model, 
we try to extend the result to Za-symmetric model. Fortunately, the known 9x9 matrix representation of the solution to YBE 
is a proper unitary operator for constructing the desired model with Z 3 symmetry. Now we give a brief introduction to the 9x9 
matrix representation of the solution to YBE which is associated with this paper. Eirstly, let us introduce the braid matrix [20] 
forw = e‘T: 



'wOOOOlOwO' 
OwOw^OOOOw^ 
OOwOwOlOO 
Ow^OwOOOOw^ 
OOlOwOwOO 
wOOOOwOlO 
OOwOlOwOO 
lOOOOwOwO 
Ow^Ow^OOOO w 


which satisfies the braid relation 


(16) 


BiBi+iBi — Bi+iBiBi+i, 


(17) 


where B{ = J®.../® B 0 /... (I is 3x3 identity matrix). 

i,i+l 

The solution Pi (0) of Yang-Baxter equation can be viewed as the parametrization of braid operators, 

Pi(0) =e'^^ = cos 0-l-i sin 0M, 

M 2 = 1 , 


( 18 ) 
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where Ti = e's (J+Bi) satisfies the Temperley-Lieb algebraic (TLA) relation[21] 

= dTi, d = \/3 

TiT.iiTi = Ti, (19) 

T,rj=TjTi, 

Then the YBE[22] means that a 3-body S-matrix can be expressed in terms of three 2-body S-matrices, i.e.: 


Ri{9i)Ri+i{62)R03) = Ri+iiWiiWi+iiSi)- 


( 20 ) 


Substituting equation (18) into equation (20), we have the constraint for three parameters 0i, 02 and 0^ : 

^ tan0i-|-tan03 

tan 02 =-i-• 

l-fi tan 01 tan 03 


( 21 ) 


When 01 = 02 = 03 = f-, the YBE turns into the braid relation with ^i(7r/3) = uiBi. Note that due to the different d in TLA 
for 4x4 and 9x9 solutions of YBE(c? = \/2 for 4 x 4), the angular relation in equation (21) for 9 x 9 is different from the angular 
relation for 4x4 in equation ( 1), 


, tan 0 ^-|-tan 03 

^ l-|-tan0'j^ tan 03 


( 22 ) 


It is well known that the physical meaning of tan 0i(or tan 0') for i = 1, 2, 3 is the velocity of a particle. The angular relation 
for 4x4 solution means the Lorentz addition of the velocity. 

By introducing the 3-cyclic SU(3) generators T® based on the principal representation of V. Kac[23] (see Supplementary), 
TLA generator can be expressed as the tensor product of nearest SU(3)-lattice sites, 

Ti = (23) 

v3 


Here [25] 


■ 0 

1 

o' 
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(24) 


Hence at each lattice site there is one SU(3) operator which can be identified with the colors blue, red and green [26]. In the 
following sections, we will make use of equation (23) to generate the topological non-trivial models and triple degeneracy. 


Ladder operators of SU(3) spin and extended Jordan-Wigner transformation. In this section, we present the ladder opera¬ 
tors of SU(3) spin and introduce the extended Jordan-Wigner transformation for SU(3) spin sites. 

For spin-2 at lattice sites expressed by SU(2) Pauli matrices, the ladder operators are 


'0 1' 


'0 o' 

0 0 

, CT = 

1 0 


Similarly, we introduce the cyclic ladder operators of SU(3) spin 



'o 1 o' 


'o 0 o' 


' 0 0 o' 

= 

0 0 0 

= 

0 0 w 

,d+ = 

0 0 0 


0 0 0 


0 0 0 



0 0 


'o 0 o' 


o 

o 

o 



'o 0 Lu' 

u~ = 

1 0 0 

,S~ = 

0 0 0 

,d- = 

0 0 0 


0 0 0 


0 0 


0 0 0 


( 25 ) 
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The above operators act on the color space to transform the three colors into each other obeying the algebraic relations 

[«+]2 = [ s +]2 = [ d +]2 = 0 , 

[n-r = [s-r = [d-r=o, 

u'^s~^ = d~, s~^d'^ = u~, d'^u'^ = s~, (26) 

s~u~ = d'^, d~s~ = ^ u~d~ = s’*", 

u'^s'^d'^+s'^d'^u'^+d'^u'^s'^ = 1. 


Unlike SU(2) spin, the SU(3) spin has 3 independent ladder operators m+, s+ and d+ to make the transition between three 
different colors, u, s and d can be expressed by it+, s+ and 

To introduce the extended Jordan-Wigner transformation for SU(3), let us review the original SU(2) J-W transformation firstly. 
J-W transformation transforms sited spin-^ operators onto spinless fermions 


at = 


.1=1 


■'n 5 


n 

,i=l 


where and am satisfy the fermionic commutation relations 

ant — dmn^ — 0- 


(27) 


(28) 


From the above transformation, it turns out that the anti commuting relations of fermions result from {erf, erf"} = (erf, CTj } = 0. 
Similarly, the extended J-W transformation for SU(3) can be defined as[24] 


where [25] 


in which it follows 


ul = 

sl = 

Dl = 


n-l 




_i=l 

n-l 




i=l 
"n—1 




,i=l 


J — 


S' — 


dt, D, = 


n-l 


nir-Fl- 

i=l 

n[T<‘>i, 


IXrj 1 


i=l 
'n-1 




i=l 


dn ) 




1 0 0 
0 w 0 
0 0 a;2 


a^(t) _ 

^3 — 


1 0 0 
0 0 
0 0 w 


XiYj=ojYjX, {{X,Y}e{U,S,Dy,i<j)- 

SiUi = dI DiSi = Ul, UiDi = Sl, 


(29) 


(30) 

(31) 


that can be checked straightforwardly. Different from the anti commuting of spinless fermions, in the exchange between the 
above operators there appear extra oj (or w^) phase factor. The physical meaning is obvious, when making exchange between 
two particles on i-th and j-th sites (i < j), the system gains an extra ui phase factor. Exchanging the two particles again, the 
system returns to the initial state. 

By introducing the linear combination of sited SU(3) operators 


F/ = uIWsI+ojDI gI = 

Fi = Ui-\-ujSi-\-uj'^Di, Gi = 

f/ = f2, g/ = Gf, f 3 = Gf = 1; 
XiYj=ujY,Xi ({X,y}e{F,G};i<j), 


the TLA generator F in equation (23) can be rewritten as 



l+wFj^Gi+i+FiG,t|_]^ 


(32) 

(33) 

(34) 

(35) 


(36) 
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Here and G} are non-local operators. Indeed, they are the generalizations of Clifford algebra which corresponds to Majorana 
fermions. Redefining 


Cli-i=F^, C2i-i=Fi; (37) 

Cl=coGl G2i=uj^Gi. (38) 


Thus the extended generator of TLA shown in equation (36) can be written in terms of the form 


T' = 






with 


aq = coQQ, (^ < j), 

a = [c^r, [Ci]' = [c/]' = 1- 


(39) 


(40) 

(41) 


For convenience we call the commutation relation shown in equation (40) w-commutation relation. This commutation relation 
can be regarded as the generalization of Majorana fermions’ anti-commuting, which is also proposed in [27]. Note that T( does 
not equal to Ti, but it also satisfies d = a/ 3 TLA in equation (19) and can be substituted into equation (18). In ID Kitaev 
model, two real Majorana operators corresponds to one complex fermion site as well as one SU(2) spin site. Similarly, the two 
w-commuting operators G 2 i-i and C 2 i coiTespond to the Lth SU(3) spin. Obviously, equation (40) looks g-commutation relation 
for = 1 in quantum algebra [28]. 


Generating Z 3 parafermionic model from YBE. From equation (18) and equation (39), we obtain the unitary solution Ri{9) 
of YBE in the form 


Ri{0) = e ‘^-|-i-sin0 


1+uj^ClCi+i+Lu^CiCl^, . (42) 

Now let us construct the Z 3 parafermionic chain based on equation (42). Substituting equation (42) into equation ( 8 ), we get 

(43) 




Similarly, we consider the nearest-neighbour interactions of Ci’s and extend equation (43) to an 2N-chain and ignore the constant 
term, the derived chain model can be expressed as: 



N N-1 

(cLlC2i+C2i-lCl^ + Y. ^2 (4C^2.+1+C2.4+,) 


(44) 


Here we emphasize that the chain possesses open boundary condition. This model is the ID Z 3 parafermionic model[24], which 
originates from the three-state Potts model[29-31]. Instead of the Z 2 parity symmetry of Kitaev model, the model in equation 
(44) possesses Z 3 symmetry. The symmetry operator is 


^=n , p'= 1 . (45) 

i=l 

Hence P is a Z 3 symmetry of the model and the eigenvalues of P is 1, w and Next we analyse the obtained model in two 
cases. 

1 . 01 > 0 , 02 = 0 . 

In this case the Hamiltonian becomes: 


Hi 


Of, ^ 

-y 010;^ E {ci-lC 2 i+C 2 i-lGl^ 
i —1 

i—1 


( 46 ) 
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TABLE I: Comparison between SU(2) and SU(3) pictures 


Type of spin sites 

SU(2) 

SU(3) 

Ladder operators 

’ <^17 ■ 

, s,"*” , . 

k ’ k ’ k 

After J-W transform. 

aj, a^. 


Quasiparticle 

72k-l=0([+ak; 

72k^i(oi[-ak). 

(Majorana fermions) 

c’L=7+‘^'7+‘^7; 

C2k=“7 +<^7+“7 • 

Commutation relation 

7i7j^-7j7i; 

cl =a,C^, (j>iy, 


7 = 1.71=7- 

[c)]3=i, [Ci]^=CiL 

Braid operator 

7=75 (i+7i7i+i); 

B, = ^uj^{l+ClCi+i+CiCl^^); 


[B.T = 1- 

[B.]S = 1. 

Cyclic operation 

S(7iS'=7i+i; 

B'7i+i-B' = -7i; 

-Bi'^(-7i)-B' = -7i+i; 
b7 (—7i+ 1) B' =7i. 

bI cl Bi=cl_^_^-, 
Blcl^^B,=uj'^Ci+iCi; 
bI (cj^Ci+iCi)B,=Cl. 


Here we note that 621-1 and 621 correspond to i-th SU(3) spin, dj = wCj; = and the 

vacuum state |0) is defined as dijO) = 0. The Hamiltonian is diagonalised and the ground state is unique. This is a trivial 
case. 


2 . 01 = 0 , 02 > 0 . 

In this case the Hamiltonian is: 


H 2 


Oh 

E ( 4C ^ 2.+1 + 62.4+ i ) 

k^l 

Of. 

i—1 


(47) 

(48) 


Here the quasiparticle at lattice can be defined as di = 621 — 02 ^ 621 + 1 . The ground states satisfy the condition di|V') = 0 
for 2 = 1,..., 1. Under the open boundary condition, it shows that the absent operatorsCi, Cj, C 2 N and C\-^ in H 2 

remain unpaired and are the symmetry operators of the Hamiltonian H 2 - Together with the w-parity operator P, C\, Cj, 
C 2 N and lead to the triple degeneracy of ground states which can be categorized according to P. The Hamiltonian 
has three degenerate ground states: \ipo), \'ipi) = d^lV’o) and \ip 2 ) = [d^]^|V’o), where rfl = C\—ujC\-f^. The three ground 
states I'lpo), I'lpi) and |i/> 2 ) possess the parity 1 , w and respectively. 


From the above discussion we see that the Z 3 parafermionic chain is natural generalization of the Z 2 Kitaev model. Next let us 
construct the topological invariant for the parafermionic chain and discuss its phase transition. In terms of Fourier transformation, 

^ m=l ^ m=l 

A, = ^ 

^ m=l ^ m=l 


m=l 


1 


1 


m=l 


(49) 
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The parafermionic operators in momentum space can be written in the following form 

N 


^A.k ~ 


1 

7^ 


m=l 


2m-1 1 


N 


cl , = ^B.k = ^ ^ e'-’‘C2rn- 


k 


y/N 




m=l 


Then the equation (44) can be expressed in momentum space as 


H = -- 


2h 


TT 


^A,k’ k 


k= —71 


^ [^A, k, ^B,k] 5 


(50) 

(51) 


(52) 


where M is 2x2 matrix. 


M = (eiuj‘^+02e-'^'j a++(eiUj+e 2 C^'^ a-- 


'0 1' 


'0 o' 

0 0 

, 0- = 

1 0 


(53) 


Here we note that [ Cg j, ] is analogous to the Majorana operators in momentum space in ID Kitaev model. If one define 
[a^, a'^, a^] as a SU(2) space, M can be regarded as a vector in XY-plane of SU(2) space with the basis and (k € [—tt, tt]), 
namely. 


M = M-a, 

= cr'*'+tT“, (T^ = —ia'^+ia~, = a'^a~—a~a'^, 

M= ^-10^+02 cos(fc+^), ^61-62 sm{k+^), 0), 



Now the topological invariant for vector M can be defined [32], 


W = 


dk 

Att 




dMft 

dk 


(54) 

(55) 

(56) 

(57) 


(58) 


Indeed, the topological invariant W means the winding number of the vector M winding around the original point in the 
first Brillouin zone. In Ref. [24], the author emphasized that the energy spectrum of parafermionic model can not be obtained 
simply by Fourier transformation due to the relation in equation (40). Here we do not expect to obtain the energy spectrum, 
but in the momentum basis of C\ and Cg j^, the topological winding number of Z 3 parafermionic chain shows the analogous 
characteristic as the Z 2 Kitaev chain. When | 02 1 > \di\, the winding number W = —1 corresponds to the topological non-trivial 
phase. When | 02 | < |0i |, the winding number VF = 0 corresponds to the topological trivial phase. In this sense, | 02 | = |^i| is the 
phase transition point. By calculating the eigenvalues of M, we can find that the “bulk gap” closes at | 02 | = |^i | where the “bulk 
gap” closes, the phase transition occurs. Thus we see from the above definition that the critical point of the phase transition 
1 ^ 2 ! = 1 ^ 1 1 coincides with the Z3 conformal field theory(CFT)[33-35]. Obviously, the above properties in our derived Z3 
parafermionic chain are very similar to ID Kitaev model. However, there are still some differences between Z2 and Z3 models. 
The critical point of Z 2 Kitaev model can be described by Ising CFT When Kitaev model is in topological phase, it appears 
Majorana zero mode with quantum dimension ^/2. While the critical point of Z3 parafermion model is also described by Z3 
parafermion CFT, but the non-abelian primary fields are not Z3 parafermion field. There are totally six different quasiparticles 
in Z3 parafermion model, three of which possess abelian fields, the vacuum I, parafermion field ip and Besides, there exist 
three types of non-abelian fields, a, a\ e, where a is the spin field and e = aip is the Fibonacci anyon with quantum qimension 
2271 . For example, the Z3 Read-Rezayi quantum Hall phase supports Fibonacci anyons, which are applicable to universal 
quantum computation[35-37]. 

Now let us discuss the generalization of the cyclic chain model from Z3 to Zat. To start with, let us introduce the irreducible 
cyclic representation of SU(N) generators. Under the N-dimensional orthonormal basis {|A) |A = 1, 2, ...TV}, akin to SU(3), the 
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ket-bra representation of 7V^ —1 SU(N) generators are 


N 


^ ^ = 12^ ...TV) 


a—1 


ry-T(ni) ry-T(n) f 771 _ 1J _ 1 ^ ^(^Tl+n l(modN)) 

i j *+i+l(moiIN) ’ 

-( 1 ) 


(59) 

(60) 


with uj root of unity uj^ = 1 and identity operator. Then the representation of Temperley-Lieb algebra Ti on (i, i+l)-th 
sites under the cyclic SU{N) operators is expressed as 


T,; = 




1 k 


. 2,2 + 1 


( 61 ) 


with X, y, m, n arbitrary certain integer given from 2 to N. Based on the algebraic relation in equation (60), it is easy to check 
that Ti in equation (61) satisfies TLA with quantum dimension d = ■\/N, 

+ = dTi, d = VN, 

TiT.iiTi = Ti, (62) 

T,Tj = TjTi, \i-j\ > 1 . 

From the view point of rational Yang-Baxterization of Temperley-Lieb algebra, when d < 2, i.e. < 4, the parameter in 
the solution of YBE is real and the corresponding +matrix is unitary and can be viewed as unitary evolution operator of a 
quantum system. That is, we can obtain Z 4 parafermion chain from YBE in an similar way as Z 3 . While d > 2, the parameter 
in the solution of YBE is imaginary(see Supplementary), hence the R is not unitary and cannot be viewed as an ideal evolution 
operator. One can still construct Zjv >4 parafermion chain, but the chain does not come from the rational Yang-Baxterization. 


Fermionic representation of Z 3 parafermionic model. In this section we express the Z 3 parafermionic chain in terms of 
fermions. Eor the 3x3 matrices in equation (25), we choose three orthonormal basis, |r) = r+ac), |p) = p+ac) and 
\b) = 6 +ac). Here |vac) represents the vacuum state, r^, and are three types of fermions and satisfy the fermionic 
commutation conditions 

{x^y} = 4y, = {a;,?/} = 0, {x,y = r,g,b), (63) 

with the constraint of the occupation number of the fermions on each site 

r^r+g^ g+b^b = 1. (64) 

Then equation (64) means the only one occupied fermion on the site for either or g'^ or Considering equation (25) the 
SU(3) ladder operators can be written as 

u'^ = r+, s’*" = uig^b, + = uP'b^r, 
u~ = pV, s~ = oJ^b^g, d~ = uir^b. 

In the basis of |r), \g) and \b), we have the relation 

w^xy^z\n) = SxyW^z\n)—w^xy\n) = SxyW^z\n}. {w,x,y,z,nG {r,g,b}) (66) 

It can be proved that the operators in equation (65) satisfy the same relations as for those in equation (26) (see Supplementary). 
In other words, the fermionic representation of the SU(3) operators also satisfy the matrix multiplication. Similarly, the nonlocal 
operators U\ and can be expressed as 


ul = 

sl = 

Dl = 


n—1 


[r/ n+ujglgi+uj'^bl bi] 

n—1 

[rln+ujgl gi+uj'^b'l k] 

"n—1 

[r/ n+ojg'lgi+oj'^b'l k] 


n-l 


Li=l 


rlgn, Un = 

Ujglbn, Sa = 
UJ^blxn, Dn = 


Yllrjn+w'^glgi+ujblbi] +r„; 


i=l 

n-l 


[r/ n+uj'^gl gi+ujb} k 

i=l 
n—1 

Ylirln+uj^ gl gi+ujblbi 


UJ‘ 


Li=l 






(67) 
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Here we do not need to require whether the commutation relation between two operators on different sites is fermionic or 
bosonic, since the non local operators U, S and D are always even power of the fermionic operators. 

Making use of the fermionic representation of W, and we find that T( given by equation (39) also satisfies T-L algebraic 

relation. Then the Z 3 parafermionic chain is rewritten as 


^ [^1 (C2tiC2. + C2i_iC'li)+02 [clC2i+l + C2,Cl^,) 

i —1 

Oh ^ 

= —^E [^1 (3ff/5i-l)+^2 (-dldi+2^ , 

i —1 


( 68 ) 


where dj = rjgi+gjbi+bjri—uj{r^_^_^gi+i+gj_^_ibi+i+bl_^^ri+i). For the topological non-trivial case 61 = 0, 02 > 0, J/ in equation 
( 68 ) shows the rotation symmetry of g^ and 5^ for each sited SU(3) spin. For the topological trivial case 9i > 0, 62 = 0, 
equation ( 68 ) shows that the ground state corresponds to the full occupation of fermion g^ on each site. Unlike the topological 
non-trivial case, there is no cyclic permutation symmetry of r\ g^ and 6 ^ on each SU(3) spin site. There are totally three 
types of parafermions on the i-th SU(3) spin site, f/, G] and uiFiGi (see equation (32)), but we only choose two of the three 
types of parafermions to represent Temperley-Lieb algebra in equation (36). Hence there are three different ways to choose the 
representation of TLA as well as the Z 3 parafermionic model. Each way corresponds to one type of ground state occupation(r^ 
or g^ or b^). 


Algebra of triple degenerate model. In previous sections, we have discussed the Z 3 parafermionic model and the physical 
consequences. There appears the triple degeneracy in ground state and the emergence of tripling corresponds to the topological 
phase of the Hamiltonian. It can be regarded as the extension of the algebra of Majorana doubling pointed out in Ref. [7] . In this 
section, we shall show the algebra of triple degeneracy at each energy level due to the 3-cyclic and give its intuitive explanation. 

Firstly let us construct 3-body Hamiltonian based on the 3-body S-matrix constrained by YBE. It is well known that the 
physical meaning of Ri{0) is 2-body S-matrix. YBE means that a 3-body S-matrix can be decomposed into three 2-body 
S-matrices in the following way 

.^123(^1, ^*2, 03) = -Ri2(0i).R23(^*2).Ri2(03) 

= i?23(03)-Rl2(02)i?23(0l). (69) 

Here we note that due to the constraint of equation (21), only two of the three parameters 0i, 62 and 63 are free. Suppose 0i and 
02 are time dependent, then the 3-body Hamiltonian can be obtained from equation ( 8 ) (see Supplementary) 

iTi23 = uj'^[a{G\G2+G,Gl)+l3{G\G3+GiGl)+-t{G\G3+G2Gl)]+Hi{G\GlGl+G,C2G3), (70) 

where a, /3, 7 and k are real parameters depending on 0i and 02 . By making inverse Jordan-Wigner transformation for SU(3) 

to transform Ci’s back into SU(3) spin sites, one can show that there are only two independent symmetry operators (see Supple¬ 
mentary) 

— L/1L/2O3L/4, 

F = a;C'JC'2C']. 

Then the complete set of the algebra for the Hamiltonian is 

[iJi23,F] = 0 , [7Ti23,r] = 0 ; (73) 

FP = wPF, p3 = 1, F^ = 1. (74) 

Here P represents the w-parity operator. Now we turn to the analysis of the degeneracy of the Hamiltonian. From equation (74), 
F transforms the common eigenstates \tjj) of iTi 23 and P to the following form: 

P\ip) = alip), (a = l,w,w^); (75) 

P[r|V:)] = w^FPlV') = aw^FlV'); (76) 

P[r2|^/))] = wF^PlV’) = awF^IV^). (77) 

Because F commutes with the Hamiltonian H 123 , the above three states have the same energy with different w-parity. As a 
consequence the Hamiltonian possesses triple degeneracy on all energy levels. In this sense, we conclude that Z 2 parity leads to 
Majorana doubling, whereas the Z3 w-parity leads to the tripling. 


(71) 

(72) 
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DISCUSSION 

Before ending the paper, we would like to make some comments and discussions. 

1) In our SU(3) models, only the three operators W, S'^ and are basic operators that are extention of the spinless fermion 

and a in SU(2) Majorana models. U, S and D can be expressed by W, and Z)^(see equation (31)). 

2) Our results can be regarded as the direct generalization of Z 2 Majorana models. The comparison between Z 2 Majorana 
fermion and Z 3 parafermion is shown in the Table I. In Z 2 case, the eigenvalues of the symmetry operators are ±1 whereas the 
eigenvalues turn into 1, ui and for Z 3 symmetry operators. We see from the Table I that the exchange of Z 3 quasiparticles Ci 
emerges an extra uj or phase factor instead of — 1 . 

3) The topological case of Z 3 parafermionic model has been obtained, which corresponds to the triple degenerate ground 
states and the non-trivial topological winding number. Although the Z 3 parafermionic model consists of SU(3) spin, the matrix 
M of the Hamiltonian in momentum space still forms SU(2)(see equation (53)). Hence we can follow the original Kitaev model 
to dehne the topological winding number for the Z 3 parafermionic model. The Z 3 parafermionic model is formed by three types 
of fermions r^, g'^ and 6 ^. For topological trivial case, the ground state corresponds to the full occupation of p-fermion. For 
topological non-trivial case, the ground state shows the cyclic rotational symmetry of r, g and &-fermions. This may be helpful 
to realize the Z 3 parafermionic model experimentally. 

4) To give an intuitive explanation about the triple degeneracy, we construct the 3-body Hamiltonian H 123 based on YBE. 
This model is the generalization of the 3-MF model pointed by Lee and Wilczek. Two independent symmetry operators F and 
P of the Hamiltonian have been found and we show that all the energy level of H123 possesses triple degeneracy. In the process 
of constructing triple degeneracy, the 3-cyclic property of the Z 3 symmetry operators plays the crucial role. 

5) Both of the two derived models are based on the 9x9 braid matrix Bi as well as the solution Ri{9) of YBE. Regarding 
Ri{9) as the time dependent unitary evolution of 2-body interaction, one can construct the local 2-body interacting Hamiltonian. 
The 2-body Hamiltonian is then extended to the desired Z 3 parafermionic chain. To obtain the 3-body Hamiltonian, we suppose 
that 3-body S-matrix can be decomposed into three 2-body S-matrices via YBE that is acceptable in low-energy physics. The 
advantage is in that the 3-body Hamiltonian inherits the w-parity symmetry from 2-body Hamiltonian. In other words, the w- 
parity symmetry of R{9) preserves the symmetry properties of Z 3 parafermionic model and 3-body Hamiltonian. This is the 
important role of YBE plays in obtaining the desired models. 

To summarize, we extend the Z 2 Kitaev model to Z 3 parafermionic model. Due to the 3-fold 3-cyclic, Z 3 model possesses 
triple degeneracy. But both Z 2 and Z 3 models have the similar topological phase transition scheme obtained from YBE. In this 
sense, YBE is a powerful tool for generating new models. How to make full use of YBE to generate more meaningful models is 
still a challenge problem. 
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I. w-CYCLIC SU(3) GENERATOR IN YANG-BAXTER EQUATION 


The w-cyclic representation of SU(3) generators are (except 


^ = 


1 0 0 
0 1 0 
0 0 1 



'o 1 o' 


'o 0 r 

rp{2) 

0 0 1 

rp(3) 

’ 1 ~ 

1 0 0 

1 0 0 

0 1 0 



1 —‘ 

o 

o 


' 0 1 o' 


'o 0 r 

II 

0 w 0 

0 0 

'l 0 o' 

II 

0 0 w 
a ;2 0 0 

'O 1 0 ' 

II 

w 0 0 

0 a ;2 0 

' 0 0 r 

II 

0 0 

0 0 (jj 

II 

0 0 a ;2 
a; 0 0 

II 

0 0 

0 w 0 


and they satisfy the following algebraic relations 

yCm) ji(n) _ ^(m-l)(j-l)2^(6) 


(w = e' 3 ), 


( 1 ) 


a = i 4 -j-l(mod 3 ), b = nn-n-l(mod3). 


The braid operator can be parametrized to yield the solution of YBE by means of Yang-Baxterization[2]. Now let us recall 
the standard method of Yang-Baxterization. The YBE reads, 


and the solution is 






G{p) = 




fiP.i') = 


tto—(i/i /2 ’ 

p+v 


/3 = 


d^-A 


( 2 ) 

(3) 

(4) 

(5) 


(2ao)^’ 

Here Ti is Temperley-Lieb algebra(TLA) generator, d represents the loop value of TLA and oq is a free parameter. If we express 
Ti in terms of the known braid operator Bi, then the solution of YBE can be obtained. In this paper, the T-L generator associated 
with braid operators Bi can be expressed as follows, 

Ti = ( 6 ) 




i,i+l • 


Ti satisfy T-L algebra with the loop value d = y/3, 

T^ = dTi,d = V3 

TiT.iiT. = Ti, 

TiTj=TjTi, 

Based on the algebraic relation in equation (1), we can verify these relations by direct calculation. 
After the replacement 


(7) 


tan 01 
tan 02 
tan 03 


P 

P 

P 


id 

y/dP^ 

id 

id 

yfdP^ 


u, 

V, {d = Vs) 


( 8 ) 

(9) 


( 10 ) 
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the solution Ri{u) is rewritten as 


with the angular relation 


Ri{0) = cos6l/+isin0(—pTi—/), 
v3 


tan 02 


tan0i+tan 63 
1+^ tan 01 tan 63 


When 01 = O 2 = 03 = the YBE turns back into the braid relation with _Ri(7r/3) = coBi. 
In comparison with d = for -R(p) related to the Bell basis [3], p = tan 0, and 


M 1 +M 3 

02 = 7 --• 

1+0103 


(11) 


( 12 ) 


II. PROOF OF OPERATOR MULTIPLICATION FOR u^,s^, IN THE EERMIONIC REPRESENTATION 

In the orthonormal basis of |r) = rljvac), \g) = (/l|vac) and \b) = 5l|vac), with the fermionic condition on each site 

{a;l,y} = 4s/, {a;l, yl} = {x, y} = 0, {x,y = r,g,b), (13) 

the operators u'^, s^,d^ are expressed as 

u'^ = r^q, s"*" = ujq^b, d~^ = ur^h^r, 

(14) 

u~ = pV, s~ = ui^b^g, d~ = Lor^b. 

Due to the constraint that the total occupation number for the fermions r, g,bis 1 on each site, 

r^r+g^ g+b^b = 1. (15) 

In the basis of |r), \g) and |6), we have the relation 

w^xy^z\n) = SxyW^z\n)—w^xy\n) = dxyw'^zln). {w,x,y,z,n € {r,g,b}) (16) 

i.e. acting two adjoint annihilation operators on the single fermionic occupation basis. Then multiplication of the operators 
u^, s^, d^ are easily checked in the basis of |r), |p) and \b), 

[u+]2 = r^gr^g = —r^r^gg = 0; 

= g^rg^r = —g^g^rr = 0 ; 

[s+]2 = to^g^bg^b = -io^g^g^bb = 0; 

[s“]^ = ojb^gb'^g = —ojb^b'^gg = 0; 

[d'^]'^ = ujb^rb^r = —ojb^b^rr = 0 ; 

[d~]‘^ = Lu'^r^br^b = —tu'^r^r^bb = 0; 

= [rlp][a;pl5] = Lor^ {1—g^ g)b = tor^b = d~] 
s'^d'^ = [ujg'^b][uP'h^r\ = g\l—b'^b)r = g^r = u~] 
d'^u'^ = [w^5lr][rlp] = rVjp = ui^b"^ g = s~; 

s~u~ = [uj'^b'^ g][g^r] = uj^b'^ {l—g^^ g)r = ui'^b^r = d'^; 
d~s~ = [ujr"^ b][uj‘^b^ g] = r'^{l—b^b)g = g = 
u~d~ = [(/lr][a;rl5] = tog^ {l—r^r)b = cog^b = s’*". 


Hence Temperley-Lieb algebraic relation represented by fermions can also be checked. 
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III. DERIVATION OF 3-BODY HAMILTONIAN 


The 3-body S-matrix constrained by YBE is 

^123(^1,^2,03) = -Ri2(0i)-R23(02)-Ri2(03) 

= i?23(03)i?12(02)i?23(0l). 

Regarding ^123 as the unitary evolution of system, one can construct the 3-body Hamiltonian 


(17) 


Hmit) 


(18) 


Here we note that due to the constraint of equation (12), there are two of the three parameters 6 * 1 , 62 and 6*3 are free Supposing 
that 6*1 and 6*2 is time dependent, we obtain 


Hm =ifi0i^4l^^r2'(0i)+ift02i?i2(0i) 


dOi 

3 


R-^\e, 


— — fi01 (■yjT'l — 1 ) — fi'02-R 12 (01) 




T 2 -I 


002 "‘■23 

i?r2'(0l) 


=-§^2 sin 6 »i sm( 6 »i+ 7 r/ 3 )w 2 (C'J'C 3 +C'iC'])+M 2 (| sin^ 6 li-|)w 2 (C'jC' 3 +C' 2 C']) 
-§^2 sin 6*1 sin( 6 li- 7 r/ 3 )(C/C' 2 ^C' 34 C'iC' 2 C' 3 )-|?i 0 iw 2 (C']'C' 2 +C'iC' 2 ^)+/ii( 6 l'i+ 02 ). 
Ignoring the constant term, we have 

a = 

P = —|/i02sin0i sin(0i+7r/3), 

7 = ft02(|sin^ 01-i), 

K = —|/i02 sin01 sin(0i—7r/3), 

for 

7Ti23 = uj'^[a{C\C 2 +C,Cl)+P{C\Ci+C,Cl) 

+7(Cl C 2 +C 2 CI )]+«;(Cf cl cl +Ci (72^3). 


(19) 


( 20 ) 


IV. SYMMETRY OPERATORS OF 7/123 


In this section, we show that there are only two independent symmetry operators of 3-body Hamiltonian H 122 , ■ Let us first 
transform equation (19) into matrix tensor product form under SU(3) Jordan-Wigner transformation 


ul = 
sl = 

Dl = 


n-l 


n[^2 i>i« 
n'pf’t 


+ 
n •> 


.1=1 


C/n = 

"> — 

dt. D,= 


n-l 


i^l 

iirt”ii 

= 1 


Li=l 


d-. 


We have 


d = , Cl = ; 

cl = ujTd®Td, C2 = Td®Td] 
cl = rd C)Td , C 2 , = ; 

cl = LoT^dT^, Ci = T^dr^- 


(21) 
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Then the Hamiltonian reads 

+7 (a;2rf^(g)r|^4wTf^(g)rp^)+K (u;T3^^^(g)Tf^+T2^^^(g)Tf^) . 

Based on the the algebraic relation in equation (1), we try to find the independent symmetry operators without linear composition. 
The general form of the operators can be expressed as 

= ( 22 ) 

Here a, b, x and y are to be determined, {a, b, x, y} € {1, 2, 3}. The condition is that each term in H 123 commutes with Ti, 

= 0 ; 

[TP(g)T^^\T^^'>(g)Tl^'>] = 0; 

= 0 ; 

= 0 . 

After direct calculation, one obtains 


x = l, a = &G {1,2,3}, t/G {1,2,3}. (23) 

Totally, there are 8 cases, 

Ti = rf^oTf^Tg = 

Tz = = Tg'^^OTg'^^ 

Tg = T2^^^(g)Tf\r7 = Tg^^^OTg^^^ 
r4 = = Tg'^^OTg^^^ 

But it is easy to check that 

Tg cx [ri]2, Te cx [Tg]", 

Tj cx [Tg]", Tg cx [r 4 ]", 
rg = r2ri,r4 = r2[ri]2. 

Then only 2 independent operators are left, 

Ti = \ (24) 

T 2 = T^^'’(E)T^^'’. (25) 

Making inverse SU(3) J-W transformation, we can redefine that 

P = r2 = CiClCsCl (26) 

r = ri = a;C'JC' 2 C']. (27) 

Here we define P as the Z 3 w-parity operator. Hence we find the two symmetry operators of Hi 23 . 
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